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O ' Abstract 

(N 

^ . We study the correlation complexity (or equivalently, the comniunication complexity) of 

, ^ I generating a bipartite quantum state p. When p is a pure state, we completely characterize the 

complexity for approximately generating p by a corresponding approximate rank, closing a gap 
left in Anibainis, Schulman, Ta-Shma, Vazirani and Wigderson (SIAM Journal on Computing, 
32(6):1570-1585, 2003). When p is a classical distribution P{x,y), we tightly characterize the 
complexity of generating P by the psd-rank, a measure recently proposed by Fiorini, Massar, 
V^ . Pokutta, Tiwary and de Wolf {STOC2012). We also present a characterization of the complexity 

^^ ' of generating a general quantum state p. 

O 



> 



*Department of Computer Science and Centre for Quantum Technologies, National University of Singapore, 2 
Science Drive 3, Singapore 117542. Email: rahul@comp.nus.edu.sg 

'Department of Electrical Engineering and Computer Science, 2260 Hayward Street, University of Michigan, Ann 
Arbor, Michigan, 48109-2121, USA. Email: shiyy@eecs.umich.edu 

^Centre for Quantum Technologies, National University of Singapore, 2 Science Drive 3, Singapore 117542. Email: 
cqtwz@nus.edu. sg 

''Department of Computer Science and Engineering and The Institute of Theoretical Computer Science and Com- 
munications, The Chinese University of Hong Kong. Email: syzhang@cse.cuhk.edu.hk 



1 Introduction 

In [5], the following basic model was studied: Two parties, called Alice and Bob, aim to generate a 
target bipartite state p € T-La ® 'Hb (Hilbert space T-La is in possession of Alice and Hilbert space 
y.B is in possession of Bob) using local quantum operations on a shared seed state a G TiA <X' Tis- 
The minimum sizq}\ of this seed state is the quantum correlation complexity of p, denoted Q{p)- 
Since Alice and Bob can always just share p itself, Q{p) is at most the number of qubits of p, so 
the correlation complexity is a sublinear complexity measure. Let {\x) \ x S [dim(?^yi)]} be the 
computational bases for TIa and let {\y) \ y G [dim('HB)]} be the computational bases for "Hb- 
We call a state p classical if its eigenvectors are the computational basis states {\x) iSi \y) \ x £ 
[dini(TiA)],y £ [dim(7^5)]}. Equivalently, it is just a classical probability distribution on the 
computational bases of T-La ®T-Lb- For a classical state p, the minimum size of a classical seed state 
is the randomized correlation complexity of p, denoted R{p). The work [5] exhibited a classical state 
p of size n with R(p) > log2(re) and Q{p) = 1. 

Above we considered the model in which Alice and Bob start with some shared state a and 
produce target state p by doing only local operations and no communication. On the other hand, 
we also consider the model in which Alice and Bob start with some tensor state a a (X" ctb and do 
some local operations and communication and produce p at the end of their protocol. The quantum 
communication complexity of p, denoted QComm(/>), is defined as the minimum number of qubits 
exchanged between Alice and Bob, such that at the end of their protocol they output p. Again, 
when p is classical, one can also define the randomized communication complexity of p, denoted 
RComm(/>), as the minimum number of bits exchanged between Alice and Bob, such that at the end 
of their protocol they output p. In [5j it is shown that for any classical state p, 

RComm(p) = R(p) = riog2rank+(P)l, 

where rank+(P) is the nonnegative ranlcl of the dim(7^^) x dim(?^B) matrix P with P{x,y) = 
((x|(8'(y|)/o(|x)(g)|y)). It turns out that for a general quantum state p it holds that QComm(/9) = Q{p) 
as well. This fact was attributed to Nayak (personal communication) in [5], and we shall see the 
reason in a later section. 

We have considered above two extreme models. Instead we can also consider the intermediate 
model where Alice and Bob start with some shared state cr and communicate between them to 
finally produce the target state p. In this case we count the size of a plus the communication as 
the resource used towards the complexity. Let us denote Q{p) to be the minimum resource used 
by any protocol which produces p. It is clearly seen that QComm(p) < Q(p) < Q(/o), and hence 
QComm(p) = Q(p) = Q(p) since QComm(/9) = Q{p). 

Our results 

In this paper, we conduct more studies on the fundamental question of bipartite state generation. 
We consider approximate versions of Q{p) defined as follows. Below F(p, p') represents the fidelity 
between p and p'. 



^The size of a quantum state a is defined to be half of the number of qubits of a. 

The nonnegative rank of a nonnegative matrix A is the smaUest number r such that A = X][=i ^i where each Ai 
is a nonnegative rank-1 matrix. 



Definition 1. Let e > 0. Let p be a quantum state in T-La ® T~Lb- Define 

Q,(p) = mm{Q(p') | F(p,p') > 1 - e; p' ^Ua® Ub}. 

Qr^(p)''=^'mm{Q(|,/.)(0|) | F(p, |<^)(<^|) > 1 - e; \cp) eUA^-HB]- 
In pp, Ambainis, Schulman, Ta-Shma, Vazirani and Wigderson showed that for any pure state 

riog2rank2e(^)l < QrW^)W) < [logs rank, (A)]. 
Above A is the diiQ{'HA) x di'ni{'HB) matrix with A{x,y) = a^^y and 

rank,(yl) ='^ min{rank(5) | \\A - B\\l < e}. 

Using Lemma m (as mentioned in the next section), one can easily construct a state {ip) G C" (8) C" 
such that rank2e(A) = 1 but ranke(^) = n/2, making the above two bounds arbitrarily far from 
each other. In this paper we show the following tight characterization. 

Theorem 1. Let e > 0. Let {\x) \ x G [dim(^^)]} be the computational bases for T-La and let 
{\y) I y G [diin{'HB)]} be the computational bases for T-Lb- Let \ijj) = '^xv'^x,y\^) *^ \y)- ^^^ ^ ^^ 
defined as A{x,y) = Ox^y Then 

Q,(|V')(Vl) = Qr^(lV')(V'l) = riog2rank2,_,2(^)l. 

Our result not only improves the bounds in pT] to optimal, but also shows that allowing a mixed 
state to approximate a pure state \^) does not help, for any \ip) and any approximation ratio e. 

Our second result is for the case of a classical state p. Previously [5] gave upper and lower 
bounds: 

- log2 rank(P) < Q(p) < min log2 rank((5). 

4 Q: QoQ=P 

Above P is given by P{x,y) = {{x\ {y\)p{\x) (S" \y)) and o is the Hadamard (i.e. entry-wise) 
product of matrices. How tight these bounds are is not clear yet, and an open question asked 
in [5] was a characterization of Q{p)- In this paper, we answer this question by showing a tight 
characterization in terms of psd-rank of P, a concept recently proposed in [3] by Fiorini, Massar, 
Pokutta, Tiwary and de Wolf. For a nonnegative matrix P, its psd-rank, denoted ra.iikpsd{P), is 
the minimum r such that there are r x r positive semi-definite matrices Cx, Dy, satisfying that 
P{x,y) = tr{CxDy). We show the following result. 

Theorem 2. Let {\x) \ x G [dim(7^^)]} be the computational bases for T-La oind let {\y) \ y G 
[dim(?^B)]} be the computational bases for T-Lb- Let 



ye[dim(Wg)] 

Let P be a [diin{7iA)] x [dim(?^B)] matrix with P{x,y) = px,y Then Q{p) = [log2 rankps(i(^)l • 



Along with the characterization R{p) = [log2 rank+(P)] (shown in [5]), it is interesting to see 
that for classical states p, randomized correlation/communication complexity is all about nonneg- 
ative rank, and the quantum correlation/communication complexity is all about the psd-rank of 
the corresponding matrix P. 

For a general quantum state p we show the following characterization of Q(/o). 

Theorem 3. Let p be a quantum state in'H.A^'HB- Let {\x) \ x £ [dim(?^^)]} be the computational 
bases for % A and let {\y) \ y G [dim(?^B)]} be the computational bases forT-Ls- ThenQ{p) = [log2 r] 
where r is the minimum number such that there exist matrices {A^ \ x £ [dim('H/i)]} and {By \ y G 
[dim(?^s)]}! each with r columns, and 

a;,a;'e[dimCM^)] 
i/,!/'G[dim{Ha)] 

The rest of the paper is organized as follows. In the next section we discuss our notation and 
some information theoretic preliminaries. In section [3l we prove Theorem [TJ In section [J] we prove 
Theorem [2] and Theorem [3l 

2 Preliminaries 

Matrix theory 

For a natural number n we let [n] represent the set {1, 2, . . . , n}. For a matrix A, we let A^ represent 
the transpose of A, A* represent the conjugate of A and A' represent the conjugate transpose of 
A. An operator A is said to be Hermitian if A^^ = A. A Hermitian operator A is said to be positive 
semi-definite if all its eigenvalues are non-negative. We will use the following fact. 

Fact 1. Let \vi), . . . , \vr) be vectors in C" for some n > 1. Then the r x r matrix A defined by 
A{i,j) = {vi\vj) is positive semi- definite. 

If A is positive semi-definite then so is A = A*. We let cri{A) > • • • > <Tn{A) denote singular 
values of A. The rank of A, denoted rank(j4), is defined to be the number of the non-zero singular 
values of A. The Frobenius norm of A is defined as ||^||2 = ^/^d^^^^)^ ^^^ ^^^ trace norm is 
defined as ||A||i = '^- ai. For e > 0, define 

rank,(A) = min{rank(S) | \\A - B\\l < e}. 

The following well-known result says that the best way to approximate A (under the Frobenius 
norm) with the least rank is by taking the large singular values part. 

Lemma 4 (Eckart- Young, [2j). Let \\A\\2 = 1 and e > 0. Then, 

k 
ranke(^) = the minimum k such that > o"j(A) > 1 — e . 

The following definition of psd-rank of a matrix was proposed in [3] . 

Definition 2 ([3j). For a matrix P G M"^™", its psd-rank, denoted rankps(i(-P); is the minimum 
number r such that there are positive semi-definite matrices Cx,Dy G C*"^^ with tr{CxDy) = 
P{x,y), \/x £ [n],y £ [m]. 



Quantum computing 

A quantum state p in Hilbert space T-i, denoted p £ Ti, is a, trace one positive semi-definite operator 
acting on Ti . The size of a state p is defined to be half the number of qubits of p. Here we take the 
factor of half because we shall talk about a correlation as a shared resource. It is consistent with the 
convention that when the two parties shares a classical correlation {X, Y), where Y = X = R for a 
r-bit random string R, we say that they share a random variable R of size r. A quantum state is 
called pure iff it is rank one. We often also identity a pure state with its unique eigenvector with non- 
zero eigenvalue. For quantum states p and a, their fidelity is defined as f{p,a) = tT{^/a^/'^pa^/'^). 

For p, lip) S 7i, we have F(/9, \ip){tp\) = y^('0|/9|V'). We define norm of IV') as |||V')|| = \/wW}- For 
a quantum state p S 'Ha^'He, we let tr-^gp represent the partial trace of p in TIa after tracing out 
T-Lb- Let p £ T-lASind\(j)) G T^yi^D^^s be such that tr-^^|^)((^| = p, then we call \(j)) a purification oi p. 
For a pure state \ip) G 'Ha'^T^b, its Schmidt decomposition is defined as \ip) = X]i=i VPi' bi) <^ \wi), 
where {\vi) £ 'Ha} are orthonormal, {\wi) G T-Lb} are orthonormal, Vi : Pi > with Y2l=iPi — 1- 
It is easily seen that r is also equal to rank(tr-^^|^)(^|) = rajik{tr-}{g\tp){ip\) and is therefore the 
same in all Schmidt decompositions of \^p). This number is also referred to as the Schmidt rank 
of {tp) and denoted S-rankd-i/;)). Sometimes we absorb the coefficients ^/pi in \vi) <^ \wi), in which 
case \vi),\wi) may not be unit vectors. The following is easily verified. 

Fact 2. Let Ua be a unitary operator on %a Oind let Ub be a unitary operator on T-Lb- Let 
\^p) e'N.Af^'HB- Then S-rank(lV')) = S-raiik((C/A Ub)\^Jj))- 

The following fact follows by considering Schmidt decomposition of the pure states involved; 
see, for example, Ex(2.81) of [1]. 

Fact 3. Let \'ip),\(j)) G T-La (^IIb be such that tr'^^|(/))((/)| = tr-^^|^)(^|. There exists a unitary 
operation U on T-Lb such that (/-Ha ®U)\tp) = \(f)), where /-^^ is the identity operator on Ha- 

The following fundamental fact is shown by Uhlmann [1]. 

Fact 4 (Uhlmann, [1]). Let p,cr £ T-La- Let |^) G T-La ® T~{.b be a purification of p and dim(^A) < 
dini(T-lB)- There exists a purification \(j)) G T-La ®'Hb of a such that F(/3,(t) = |((/)|'(/')|. 

We define the approximate Schmidt rank as follows. 

Definition 3. Let e > 0. Let |^) be a pure state in T-La ®T-Lb- Define 

S-rank.dV')) =^ min{S-rank(|,/.)) | |0) G ^a ® ^B and F(|V^)(V^|, |0)(0|) > 1 - e}. 

Let {\x) I X G [dim(?^A)]} be the computational bases for "H^ and let {\y) \ y G [dim('?^s)]} be the 
computational bases for T-Lb- We define linear map vecinv which takes vectors in T-La'^'Hb and maps 
them to operators from TIb to "Ha- For all x G [dim.{TiA)],y G [diui{TiB)] define vecinv(|x) (8) \y)) = 
\x){y\ and extend to all vectors in T-La (X" T-Lb by linearity. For \ip) G T-La ^SD T-Lb, it is easily seen that 
|||(/))|| = ||vecinv(|^))||2. 

In the following sections we assume Hilbert spaces T-La,T~{.Ai,T~{.a2 etc. are possessed by Alice 
and Hilbert spaces T-Lb-,T-Lbi-,T-Lb2 ^tc. are possessed by Bob. 

We start by showing the following key lemma which we will use many times in the following 
sections. 



Lemma 5. Let p he a quantum state in T-La ® T~Lb- Then, 

Q(p)=^min {[log2(s-rank(|V')))l | p = tr«^ ^^WsjV'XV'l}- 



def 



Proof. Let r = min^^^ ,^g J [ loga (s-rank(lV')))] | p = tr^^^^iKfljV')(V'|}- We first show Q(p) < 
r. Let \ip) be such that 

r = [log2 (S-rankdV-)))] and p = trn^^^nB^ 
Let t = S-rank( I '(/')). Let {ip) have a Schmidt decomposition 



E 



fi CjO ■Wi 



Let Alice and Bob start with the state 



and transform |(/>) to \^) using local unitary transformations. This shows that Q(p) < [log2 i] = r. 

For the other direction let s = Q{p)- Let Alice and Bob start with the seed state a and 
apply local completely positive trace preserving maps $a, ^b respectively to produce p. Let us 

def 

assume without loss of generality that the number of qubits of a a = tr-^^o" is at most s. Let 
a A = Yli=i '3^j|^i)(^j|i where Oj > is the i-th eigenvalue of a a with eigenvector \vi). Define 



2' 
def 



^\/ai • \vi) (g) \vi 



i=l 

and let 

2= 

!<;/'') = X]\/«i • \vi) <X> \wi) 

i=l 

be a purification of cr, where Vz : |fi) G Ha and |tt;j) G "Hb T-Ib2- 

Now consider the following operations by Alice and Bob. They start with the shared state \(j)). 
Bob using local unitary (after attaching ancilla |0) if needed) transforms \(j)) to \(j)') (Bob can do 
this follows from Fact [3]). Alice and Bob then simulate their maps ^a,^b on a by local unitaries 
(each after attaching ancilla |0) if needed on their parts; such a simulation is a standard fact, please 
refer to [4]) and finally produce a purification \6) E TIa (8) V-Ai <S5 "Hb (8) V-Bi of p. Since Alice and 
Bob, using local unitary operations and attaching ancilla |0), transform |0) to \6), we have (using 
FactE]) 2^* > S-rank(|(/))) = S-raiik(|6l)). This shows that r < s. D 

The following lemma is credited to Nayak (personal communication) in [5]. 

Lemma 6. For a quantum state p £ T-La x ^_b , Qip) = QComm(p). 



Proof. Clearly Q{p) > QComm(p). For the other direction let r = QComm(p). Let Alice and Bob 
start with the state aA® ctb G T-La (8 T~{-B, do local quantum operations, communicate r qubits and 
at the end output p. This protocol can be converted into another protocol where Alice and Bob 
start with a purification {(p) G T-Ia iX" T-Iai ® T~{-B '^ T~{-Bi of aA C3 ctb (with S-rank(|(/))) = 1), do local 
unitaries, exchange r qubits and at the end output a purification \ip) £ 'Ha'^'Hax 'iS'Tis ®'Hbx of p. 
Since local unitaries do not increase the Schmidt rank of the shared state and exchanging r qubits 
increases the Schmidt rank by a factor at most 2'' (since the rank of the marginal state possessed 
by Alice increases by at most a factor 2 on receiving a qubit from Bob, and similarly for Bob on 
receiving a qubit from Alice), we have S-rankd-;/;)) < 2'". Hence from LemmaO Q(/j) <r. D 

3 Correlation complexity of approximating a pure state 

In this section we prove Theorem [H We start by first characterizing the approximate Schmidt rank. 

Lemma 7. Let e > 0. Let \tjj) he a pure state in T-La ® T~Lb with a Schmidt decomposition 1-0) = 
Z]i=i \/Pi ■ \'^i) ® l^j) (with pi > P2 > ■ ■ ■ > Pr > o,nd Y^\^iPi = Ij- Let r' he the minimum 
number such that X][=iP« ^ (1 ~ c)^- Then r' = S-rankgd-^)). 

Proof. We will first show that r' > S-rankedV'))- Let q = '^l^iPi- Define 

1 '■' 
\(p) = —p-y^ VPi ■ \Vi) O \Wi). 

Then F{\'ip){ip\,\(j)){(j)\) = \{tp\(t))\ = ^ > 1 - e. Clearly S-rank(|(/))) = r' and hence r' > 
S-rank.dV^)). 

Now we will show r' < S-ranke(lV'))- Let s = S-rank^dV'))- Let 16*) G TIa ^T-Ib be a pure state 
such that \{6\ij)\ = f{\Tp){Tp\,\0){0\) > 1 - e and S-rank(|6l)) = s. Without loss of generality (by 

multiplying appropriate phase to \9)) let us assume that /3 = {(p\0) is real. Let 

s 

be a Schmidt decomposition of \6). Define 

r s 

del V^ ^ , , , , , „ def 



"^ = '^VPi-\'"i){'^i\ and B = f3 ■'^^i-\vi){wi\. 

i=l 1=1 

Note that A = vecinvd-^)) and B = vecinv(|0')). Since {vi} and {wi} are orthonormal, {^/Pi} form 
the singular values of A. Similarly {/3 • ^/qi} form the singular values of B. Now, 

l-{l-ef>l-(3^ 

= \\\4>)f + \\\9')f-2{9'\<P) = m-\9')\\' 

= ||vecinv(|0) - \9'))\\l = ||vecinv(|0)) - vecinv(|0))||i 

= \\A-B\\l 



Hence from Lemma IH S-rank(^) = rank(i3) > r'. 



D 



We can now get the desired characterization for Q^''^'^(|'0)('0|). 

Theorem 8. . Let e > 0. Let |^) be a pure state in T-La ® T~iB with a Schmidt decomposition 
IV') = Ya=i \fVi ■ bi) "^ \wi) (with pi> P2> ■■■ >Pr > and Yll=iPi = '^)- ^ei A = Ya=i \fVi ■ 
\vi){w,\=yeany{\ip)). Then, QP^'%\^p) {ip\) = [logs raiik2,_,2 (^)l . 

Proof. From the definitions and Leninia[5]it is clear that Q'^^^'^ {\^p) {tl^\) = [log2 S-rank^d^))]. Also 
from Lemma m and Lemma [7| it follows that S-rank^d';/')) = rank2^_^2(^) (by noting that {y/pi} 
form singular values of ^) . D 

The following lemma shows a monotonicity property for the approximate Schmidt rank. 

Lemma 9. Let e > 0. Let {ip) be a pure state in T-La "X" ^b- Let \6) be a pure state in T-Lai 'S> T-Lbi ■ 
Then, 

S-rank.dV') ® \9)) > S-rank.dV')). 

Hence from Lemma 

Proof Let |^) = EI=i V^" kl> ® 1^}) (with pi > • • • > p^ > 0) and \e) = El=i ^/^- luf) <S) \vf) be 
some Schmidt decompositions of {ip) and \0) respectively. Then 



*j 



def 



Fix a minimal set S C [r] x [s] with ^/^ i^^sPilJ — ^ ~ ^- Let r' = S-rank^dV'))- We will show 
\S\ > r'. Assume for contradiction |5| < r' — 1. Let Si = {i \ 3j such that (i, j) G S}, then 
\Si\ < \S\ <r' -1. We have 

^ PiQj < ^ Pi < Pi H h P|Si| < 1 - e, 

{i,j)&S ieSi 

where the first inequality is because Ylj-u j)£s1j — ^ ^^^ ^^^ ^' ^^^ second inequality is because 
Pi's are in the non-increasing order, and the last one is by the definition of S-rank^dip)) = r', the 
smallest number such that Pi + ■ ■ ■ + Pr' > 1 — e (from Lemma [Tj) • This contradicts the way we 
picked S and hence 

S-rank.dV') 16*)) = \S\ >r' = S-rank.dV')). D 

Theorem 10. Let e > 0. Let |V) be a pure state in Ua^Ub- Then, Qe(|V')(V'l) = Qr''''(IV')(V'l)- 

Proof. By definition, we have Qe(|V')(''Al) ^ Qe"'~'^(lV')(''Al)- Now consider the other direction. By 
the definition of Q^{\'il)) {ip\) , there exists a p G T-La ® Hb such that 

Q,(|^)(^|) = Q{p) and F(p, |^)(^|) > 1 - e. (1) 

By Lemma O there exists a purification \(p) in Ha iX" Hai (X" Hb (X" Hbi of p with 

Q{p) = riog2 S-raiik(|,/.)(,/.|)l = Q(|0)(0|). (2) 

Without loss of generality, we can assume that dim('?^^^ (g) Hbi) > dim(?^^ (g) Hb) (otherwise 
we can attach |0) to \(j)) appropriately). Now by Uhlmann's Theorem, there exists a purification 

7 



IV'') £ T~iA ® 'Hai (X" T-Ib fS) T-Ibi of \tp){tp\ such that |((/)|'(/'')l ^ 1 — £• Since \ip) is a pure state, 
1^') = 1^) O 16*) for some |6l) G ?^Ai «> ^Bi- Therefore, 

Qr''^(|^)(^|) < Qr'''^(|^')(V''l) (from LemmaED 

<Qi\(t>){(t>\) (from the definition of QP''''%\i;'){ij'\)) 

= Qip) (from Eq. i^) 

= Q.(|V')(V'I) • (fromEq. ©) 

D 

Theorem [1] now follows immediately by combining Theorem [8] and Theorem [10] and noting that 
the matrix A as defined in the statement of Theorem [T] is yec\nv{\^)). 

4 Correlation complexity of a quantum state 

In this section we show characterizations of correlation complexities for general quantum states and 
also for classical states and prove Theorem [2] and Theorem [3l We start with the following lemma. 

Lemma 11. Let p be a quantum state in T-La®T~Lb- Let {\x) \ x E [dim(?{A)]} he the computational 
bases for T-La o,nd let {\y) \ y G [dim(^B)]} be the computational bases for T-Lb- There exists a 
purification \tjj) of p, with S-raiik(|'i/')) = r, if and only if there exist matrices {A^ \ x G [dim('?^/i)]} 
and {By \ y € [dim(^B)]}, each with r columns, such that 

p= Yl \x){x'\^\y){y'\-tr[{AlA,f{Bl,By)) . 

a;,a;'6[dim{■Hy^)l 
i/,!/'6[dimCHB)l 

Proof. We first show the 'only if ' implication. Let \i/j) he a purification of pin T-La'S'T-Lai'^'H b'^T'I-Bi- 
Let S-raiik(l^)) = r. Consider a Schmidt decomposition of \ip). 

IV') = ^ b'> ® k*> = ^ Y. k)®l4)p Y. \y)'^H) 

«=1 i=l \x£[dim{HA)] / \s/e[dim(Wfl)] 

Above for any i,x,y, the vectors |w*), \w^), |w^), \Wy) are not necessarily unit vectors. Consider 

= tr^A,®wsi [y(y\x)0K)] ^lY\y)^\wi)]] 5^f J^(x'i®(f^,ij iY^y'\^H'\ 

= ^^nA,^nB,Y\ Yl \x)i^'\^K)ivi'\]^\ Yl \y)iy'\^H)H'\ 

i,j=l \x,x'e[dim{HA)] / \y,y'&ldim{HB)] 

= E( E \x){x'\-tr{\vi){vi,\)\0[ Y \y){y'\-tri\wl){wi,\) 

i,j=l \x,x'£[dim(HA)] / \y,S/'e[dim(«s)] 



Y i^)^^'i ® \y^^y'\ Yl ^^^'I4> • (<ky) 



X ,x' eldimCH j^)] \i,j = l 

V,y'eldim('Hg)] 



For each x G [dim(?^/i)], let us define matrices A^ = (|vj.), |f^), . . . , |w![.)). Similarly for each 
y G [d\u\iJ-LB)]-, let us define matrices By = {\Wy), \Wy), . . . , \Wy)). Then from above, 



E l^)(^'l ® \y)(y'\ ■ t^ (i4'^-fiBl'^y^ 



3:,a:'e[dim(?i^)l 

y.v'eldim(ng)] 



Next we show the 'if implication. Let there exist matrices {A^ \ x G [dim('H^)]} and {By \ y G 
[dim('H_B)]}, each with r columns, such that 

P= E \^){x'\^\y){y'\-tT[{AlA^f{Bl,By)) . 

y.y'e\dim(HB)] 

For i G [r], let |f^) be the i-th column of Ax and let \Wy) be the i-th column of By. Define 



def 

i=l \ X / \ y 



It is clear that S-rankd?/))) = r. We can check, by analogous calculations as above, that 

By combining Lemma[5]and Lemma [TT] we immediately get Theorem[3l We now show Theorem[2] 
which we restate below for convenience. 

Theorem 12. Let {\x) \ x G [dim('HA)]} ^e the computational bases for %a <ind let {\y) \ y G 
[dim(7^B)]} be the computational bases for T-Lb- Let 



P 



XI p^^y l^)^^l ®\y)iy\ ■ 



a:G[dimCH4)l 
!/e[dim{Ha)] 



Let P be a [dim(?^^)] x [dim(?^s)] matrix with P{x,y) = Px,y Then Q{p) = [log2 rankpsd(-P)l • 

Proof. We will first show Q(p) < [log2 rankps(i(^)l • Let r = za.iikpsd{P)- We will exhibit a 
purification |^) of p with S-rank(|^)) = r. This combined with Lemma[5]will show Q{p) < [log2 r~\ . 
Let Cx, Dy G C''^'^ be positive semi-definite matrices with tr{CxDy) = P{x, y), Vx G [dim(7^^)], y G 
[dim^g]. For i G [r], let \Vx) be the i-th column of \/C^ and let \Wy) be the i-th column of \jL)y. 
Define |^) in 'Ka ® ^A ® 'Ha^ ® ^s ® 'Hb ^ 'Hbi as follows. 

i=l \xe[dim{HA)] J \s/e[dim(Ws)] 

It is clear that S-rank(|^)) = r. Also, 

ye[dim«s] 

= E |x)(x|®|y)(y| •tr(C^L»y) = p. 



Note that Alice and Bob after sharing |'0) can either just output their first registers or measure 
their first registers in their respective computational bases to obtain p. 

Now we will show Q{p) > [log2 rankps(i(^)l • Let IV') G T-La ® T~{-Ai (^ ^s <^ ^^Bl be a purification 
of p with S-rank(l^)) = r and Q{p) = [log2 r] , as guaranteed by Lemma O We will show r > 
rankpsd(-P) and this will show the desired. Let 

i=l \xe[dim{HA)] / \j/G[dimCHi3)] 

For all X € [dim 7^^], define rxr matrices Cx such that Cx{j, i) = {vi\vl.) for all i,j G [r]. Similarly 
for all y G [dim(?^s)]) define rxr matrices Dy such that Dy{i,j) = {wy\Wy) for all i,j G [r]. From 
Factdl Cx,Dy are positive semi-definite for all x G [dim('?^yi)] and for all y G [dim(^B)]. Consider 

= Y. I^)(^l®ly)(yl ( Z]^^^l^^)-Kl^y) 

a:S[dim(K^)] \«J = 1 

ye[dim(«s)] 

= Y^ \x){x\®\y){y\ -triCxDy) . 

ye[dim(Hs)] 

Therefore for all x G [dim(^^)] and for all y G [dim('WB)] we have px,y = P{x,y) = tr{CxDy). 
Hence zankpsaiP) ^r. D 
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